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Abstract. - An adiabatic piston, separating two granular gases prepared in the same macro- 
scopic state, is found to eventually collapse to one of the sides. This new instability is explained 
by a simple macroscopic theory which is furthermore in qualitative agreement with hard disk 
molecular dynamics. 



The problem of the so-called adiabatic piston was described in 1960 by Callen [1] and has 
received a lot of attention recently [2, 3]. This construction consists of a cylinder containing 
two gases separated by an adiabatic piston. It is of interest because the principle of maximum 
entropy does not predict the equilibrium state [1, 2]. In particular any state with gases left and 
right at equilibrium with the same pressure is a possible equilibrium state of the compound 
system, even if the gases are at different temperatures. The degeneracy is however lifted when 
one moves away from the macroscopic limit. When fluctuations are taken into account, the 
set-up becomes a Brownian motor [4] and the piston moves towards the cold chamber until the 
system relaxes to "full" equilibrium with equal temperatures and pressures in both gases. In 
this letter, we raise a different question: what happens when the gases are granular, i.e., they 
gradually loose energy due to dissipative collisions. Our motivation is twofold. On the one 
hand, granular matter has recently been the object of intensive theoretical and experimental 
research [5]. Somewhat surprisingly, freely moving boundaries, which can be easily realized 
in experiment, have not received much attention. On the other hand, we will show that the 
situation is in a sense more dramatic than in the above mentioned non-dissipative case. Of 
particular interest to us is the case when both gases are prepared in the same macroscopic 
state. One expects that the piston will not move because of the left-right symmetry, which is 
preserved as the gases cool on both sides by dissipative collisions. However, our analysis will 
reveal that this state is unstable. Fluctuations or small disturbances will induce a motion of 
the piston which amplifies and finally leads to a full collapse of one of the gases. The origin 
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Fig. 1. - Time evolution of the system for the case of Ni = iV 2 = 100, Ti(0) = T 2 (0) = 1, a = 
0.99, a = 1 and M = 2000m, for t = 0, t = 10 6 and t = 10 7 . Piston dimensions are given in units of 
a. Collapse of the piston to the left side is observed on the rightmost panel. 



of the instability is simple: the motion of the piston in a given direction will compress the 
gas, increasing the energy dissipated by collisions and hence cool the gas, leading to a further 
decrease of the pressure exerted from this side. This simple intuition is confirmed by extensive 
molecular dynamics and a stability analysis of macroscopic equations of motion. 

To illustrate the type of behavior that is observed, we start by presenting some results from 
molecular dynamics simulations for a two-dimensional system. We consider a rectangular 
container of size L x x L y , separated in two compartments by a piston of vertical length L y and 
mass M. The latter can slide without friction along the horizontal x-direction. Its position will 
be denoted by x e [0, L x ). The compartments left and right, whose corresponding variables 
will be identified by a subscript i = 1,2, respectively, contain granular hard disk gases of N± 
and N 2 particles, of diameter a and mass m. The disks move freely between collisions, while 
upon collision, total momentum is conserved but a fraction 1 — a 2 of their kinetic energy is 
dissipated, where a is the so-called coefficient of normal restitution. The piston itself does not 
conduct heat and undergoes perfectly elastic collisions with the gas particles. 

Figure 1 shows the typical behavior, starting from a symmetric situation with the piston in 
the middle, x — L x /2, and equal densities and temperatures in left and right compartment. 
The parameter values correspond to small dissipation, a = 0.99, and initial low density, 
m = ri2 = 0.04. Clearly, the initial symmetry is somehow broken leading to a final collapse 
of the piston to one side. More details of the dynamics are shown in fig. 2. a, including the 
piston position and temperature in both reservoirs versus time. The black dots correspond 
to the configurations shown in fig. 1. At very short times, the piston performs small-scale 
random deviations around the starting position, as a result of the fluctuations in the pressure 
(fig. 2. a, inset). As these deviations increase, the motion develops into more regular oscillations 
combined with a net motion towards one of the sides (towards the left chamber in fig.l). The 
piston approaches one of the side walls while the gases cool down, and the oscillations shift to 
lower frequency and amplitude. In the final stage the piston completely collapses, in a finite 
time, to one of the sides, with the particles reaching a close-packing structure, as seen in the 
right panel of fig. 1. From then on, the piston stops moving. In fig. 2.b, we reproduce the results 
for higher dissipation, a = 0.85, with as most notable difference the complete suppression of 
the oscillations. Note that the direction in which the piston moves is completely determined by 
initial fluctuations. This is confirmed by the fact that the collapse is, in different realizations, 
equally often to each side, while the characteristic "onset" time for the observed instability 
also changes from one realization to the other. 

To give a theoretical explanation of the observed behavior, we proceed with the derivation 
and analysis of macroscopic equations of motion. The energy in each chamber changes due 
to two mechanisms. The first one is compression or expansion of the gas, and it is given by 
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Fig. 2. - Piston position and temperatures of the chambers versus time. Left column: parameter set 
as in fig. 1. Marked with black dots are the configurations plotted in fig. f . Right column, the same 
parameter set except a — 0.85. 



the expression dUi = —PidVi. The second mechanism is the cooling due to inelastic collisions 
that is described by Haff's law [6], i.e., the amount of energy dUi/dt\c dissipated per unit 
time is proportional to the number of collisions multiplied by the average energy lost per 
collision. More precisely at low densities it has the form dUi/dt\c = A(l — a 2 )nikBNiT^ 2 . 

Here n, = Ni/Vi is the number density, and A = 7T 1 / 2 k^ 2 'amT 1 ! 2 a constant. The volumes 
are V\ = xL y and V 2 — {L x — x)L y . We furthermore employ the equation Ui — NiksTi, exact 
for hard disks, and assume, for simplicity, the ideal gas law for the pressure, Pi — NiksTi/Vi. 
Therefore, the variation of the temperature in each chamber is written as: 

<ffi = T lV A(l-a 2 )N l T3/2 

dt x L y x 1 

m = T 2 v ^(l-a 2 )iV2 r 3/2 (1) 

dt L x - x L y (L x - x) 2 

where v is the speed of the piston, v = dx/dt. Equations (1) lead after simplification to the 
following set of evolution equations for the temperatures: 

-7i2? /2 = -72T* /2 = f t (L x -x)T 2 , (2) 

where the constants -fi = A(l — a 2 )Ni/L y have been introduced. This set of equations is 
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closed by the equation of motion for the piston: 

M^=L y (P 1 -P 2 ) = N 1 k B ^-N 2 k B -^. (3) 

(XL X X 



An analytic treatment of the above set of nonlinear equations appears to be difficult. 
However, the question of stability of a (time-dependent) reference state can be studied by 
analysis of the linearized equations. In particular the (in)stability of the symmetric state 
discussed in the above simulations can be investigated. More precisely, considering N\ = 
N 2 = N (hence 71 = 72 = 7), we note that x{t) = x(t = 0) = L x /2 and T x {t) = T 2 (t) 
T (1 + jtV%/L x y = T(t), where To is the common initial temperature, is an exact time 
dependent solution of eqs. (2)-(3) describing the symmetric cooling of the system with fixed 
piston position. To study the stability of this state, we substitute x(t) = L x /2 + £(t) and 
Ti(t) = T{t) + §i{t) into these equations, and linearize with respect to £(i) and $i(i) and their 
derivatives. We thus find: 

-\int)^m = ±±T (tm + \L,f, ( 4) 
§k§ = -i T(tm +*(«>-««). 

where the plus sign is for i = 1 and minus sign for i — 2. Integration of (4) leads to the 
following explicit solution: 

Mt) = -Mt) = ~T(tm + ^T(tf/ 2 fdt'ttt'). (5) 

L>x L x Jo 

corresponding to the relevant initial conditions $i(0) = ^2(0) =0. Substitution of these results 
in the eq. (4) gives an integro-diffcrcntial equation for £(f). By differentiating this equation 
with respect to time, and switching to a new time variable s = 1 + "f^/Tot/Lx one finds after 
some trivial rearrangements the following third order differential equation: 

where T = Mj 2 /(l6Nk B ) = tt(1 - a 2 ) 2 nL x Ma 2 /(32L y m) and n = 2N/{L x L y ). The new 
time variable s is proportional to t, and is not the time measured in units of particle collisions 
as it is usually done in studies of granular fluids [5]. Equation (6) is homogeneous in s, and 
therefore its general solution is given by: £(s) = + a 2 s^ 2 + a 3 s M3 , where the coefficients 

a, are determined by the initial conditions (initial position, speed and acceleration) and where 
the ,u's arc the roots of the equation 2T/J, 3 + 2(1 - T)fj, -1 = 0. For < T < 1, this third 
order polynomial in \i has one positive real and two complex conjugated roots. The positivity 
of the real root indicates the announced instability of the reference state, while the presence 
of complex conjugate roots show that the instability is accompanied by oscillations. 

In the simulations, T is actually very small: for instance, with a = 0.99, n — 0.04, 
M = 2000m and L y = 50, L x = 100 (parameter set in fig. 1), one finds V ~ 6.2 x 10~ 3 . 
So for comparison, we can use a small T expansion for the roots: /ii ~ ^[1 + |T] and 
[i 2 = ^3 ~ — j[l + |r] + -j=[l — j|r]. As initial conditions (remembering that t — is 
equivalent with s — 1), we consider the piston exactly in the middle, £(s = 1) = 0, with 
initial speed v, d£/ds(s = 1) = L x v / '(7\/7o), and zero acceleration, d 2 t;/ds 2 (s = 1) = 0. After 
evaluation of the coefficients a\,a 2 ,as, returning to the original time variable, and introducing 
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Fig. 3. - Left: Piston position for the parameter set of fig. 1 as a function of ln(£) revealing oscillations 
in agreement with eq. (7). Right: Comparison of the simulation results with the analytic solution of 
cqs. (1) and (3), with initial condition corresponding to the dot marked in black. 



the time to = L x /{^^/Tq) and the phase tan</> ~ 2/(3-\/r), O<0<7r/2, we finally obtain: 



3 A(l -a 2 ) Mv L x 
32 V B ~^/%T~ y 



(i + tAo) 1/2 



^=(1 + t/to)- 1 / 4 cos{0 + -L l n (l + t/t )} 



(7) 



This expression for £, although derived in the linear approximation, reproduces many of the 
features of the simulation results for small dissipation: any initial speed, however small, 
results in a superposition of a systematic motion with oscillations of decreasing amplitude 
and increasing wavelength. Quantitative agreement between theory and simulations can be 
achieved in the following way. First, according to the above result, eq. (7), the oscillations of 
the piston position plotted as a function of ln(i) have, for t large enough, a constant frequency, 
proportional to 1/y/T. Performing this plot from the numerical data, see fig. 3, we find 
r = 5.9 x 10~ 3 , while the analytical value is T = 6.2 x 10~ 3 . From the same plot we can also 
extract the value of 7 by equating the difference between the two maxima to 27r. Here we 
obtain differences of the order of 3 to 9% between analytic and numerical values. Secondly, 
while we cannot fit the solution (7) directly to the simulations starting at time t = 0, because 
the initial phase there is dominated by the fluctuations, we can consider an initial condition 
at a time t after the random motion of the piston. In fig. 3, right panel, we have included 
the analytic solution eq. (7) for the initial condition marked as dot, and observe very good 
agreement with the simulations. 

The above linear theory is only valid for not too large times: the deviation £ must be 
small, which means that the approximation is valid in the time regime vt -C WOav/v where 
v = (2/cbTq/to) 1 / 2 . In particular it cannot describe the final relaxation of the piston. We 
now proceed to show that the set of differential equations (2) and (3) predicts a final collapse 
of the piston. Let us assume that the piston is already far displaced to one of the sides, say 
the left compartment, so that x -C L x . Furthermore, we assume that the temperature in the 
right compartment does not change very much during the time interval in which the collapse 
is taking place. Hence, the force exerted by the right compartment remains approximately 
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Fig. 4. - Left: Solution of eqs. (2) and (3) where the pressure of chamber 2 is replaced by a constant 
value. The solid line corresponds to the solution of eq. (8), while the dotted line corresponds to a 
solution which initial speed does not match eq. (8). Right: simulation results when the right chamber 
is replaced by a large one that provides an approximately constant pressure. 



constant. With the replacement F ~ P 2 L y a constant in (3), one verifies by inspection that: 

r 1 (t)-r 1 (o)^i--j , x{t) - -q^W), t c = _ (8) 

is an exact solution of (2) and (3). In words, both the piston position and the temperature 
in the collapsing compartment approach their final values x — and T — in a finite time 
t c following a parabolic trajectory. This solution, eq. (8), implies specific values for initial 
position and speed of the piston. The typical solution, which we obtained by a numerical 
solution assuming a fixed force F but with an initial speed that does not match eq. (8), 
displays the same basic features but with superimposed oscillations, see fig. 4. a. Furthermore 
the frequency of these oscillations diverges as one approaches the final collapse. Linearization 
around the solution (8) predicts a divergence as l/^/x. A quantitative comparison with 
simulations is difficult because the physical assumptions behind eqs. (2) and (3) break down 
close to collapse. Nevertheless the basic features, collapse in a finite time accompanied by 
oscillations of increasing frequency are clearly observed, see fig. 4.b. 

We have focused the above discussion on the surprising properties of the symmetric case. We 
however also expect some counter-intuitive behavior in the non-symmetric case. For example, 
consider as initial condition the piston in the middle with equal temperatures left and right 
(Ti(0) = T 2 (0) = 1), but higher density (hence pressure) left (m = 0.02 > n 2 = 0.008). The 
time evolution of the piston position is shown in fig. 5. The initial displacement of the piston 
is as expected to the right side, i.e., towards the compartment with low pressure. However, 
after some time the piston reverses its motion to end, after a series of oscillations, collapsing 
to the left side. The surprising collapse of the high pressure compartment appears to occur 
whenever the initial density difference is large enough. For the "complementary" case, equal 
density (rii — n 2 ), but higher temperature (hence pressure) left (Ti(0) > T 2 (0)), one observes 
"normal" behavior with the piston eventually collapsing to the low pressure side. 

To conclude, the macroscopic equations (2)-(3), based on simple physical considerations, 
capture surprisingly well the qualitative and quantitative features observed in the simulations. 
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Fig. 5. - Piston motion for the case Ni = 50, N 2 = 20, Ti(0) = T 2 (0) = 1, L x = 100, L v = 50, M = 
2000m, a = 0.95 resulting that -Pi(O) > -p2(0), but eventually the piston moves towards the left 
chamber. Solid line is the simulation result and dotted line is the solution of the macroscopic equations. 

Agreement between theory and simulations improves significantly when one includes finite 
size corrections [7] both for collision frequency in Haff's law and for the equation of state. 
Other improvements will clearly require a much more involved theory taking into account, 
e.g., the appearance of hydrodynamic flow and inhomogeneities in density and temperature 
in the compartments, dissipativc effects at the boundary between piston and particles [8], 
deviations from Haff's law [9], and thermal fluctuations. On the experimental side, we note 
that there is some similarity between the phenomena observed here with the accumulation of 
granular gas initially distributed over several compartments, into a single one [10]. 
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